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According to classical Boltzmannian thermodynamics, the efficiency of a cyclic machine is strictly 
lower than one. Such a result is a straightforward consequence of the second principle of thermo- 
dynamics. Recent advances in the study of the thermodynamics of long-range interacting system 
report however on a rather intricate zoology of peculiar behaviors, which are occasionally in contrast 
with customarily accepted scenarios, dueling with intuition and common sense. In this paper, a ther- 
modynamical cycle is assembled for an ideal device working with non-Boltzmaniann long-range fluid 
and operating in contact with two thermal reservoirs. The system is analytically shown to violate 
the second principle of thermodynamics, a phenomenon that ultimately relates to the existence of 
regions with negative kinetic specific heat, in the canonical ensemble for the system under scrutiny. 
We argue that the validity of the second principle of thermodynamics can be possibly restored, by 
revisiting the definition of canonical ensemble, as well as the Fourier law of heat transport, and 
consequently relaxing the constraint on the maximal efficiency as imposed by the Carnot theorem. 

PACS numbers: 



The second principle of thermodynamics represents a 
fundamental milestone of physics It is a principle, 
and as such cannot be proven on a rigorous mathemat- 
ical basis. Its validity is however supported by an im- 
pressively large gallery of experimental facts and con- 
sequently assurged to the status of inviolable universal 
law. It admits different formulations, the first of which 
is credited to the German physicist Rudolf Clausius Q: 
No process is possible whose sole result is the transfer of 
heat from a body of lower temperature to a body of higher 
temperature. This in turn traces back to the concept 
of irreversibility: the heat cannot flow from cold to hot 
regions without external work being performed on the 
system, as e.g. it happens in a refrigerator. 

An alternative formulation of the second principle due 
to the Irish physicist Lord Kelvin Q recites instead: No 
process is possible in which the sole result is the absorp- 
tion of heat from a reservoir and its complete conversion 
into work. In practice, one cannot entirely transform the 
heat absorbed by a reservoir into work, without dissi- 
pating part of it. Rephrasing the latter statement: it 
is thermodynamically impossible to obtain a 100 % effi- 
ciency from a cyclic heat engine. 

Indeed, Carnot's theorem Q limits the maximum ef- 
ficiency rj for any possible engine, and can be seen as a 
direct byproduct of the second principle. All irreversible 
heat engines operating between two heat reservoirs, re- 
spectively characterized by thermodynamic temperatures 
Tc (cold) and Tjj (hot), are less efficient than a Carnot 
engine working in contact with the very same reservoirs. 
The Carnot maximal efficiency can be straightforwardly 
quantified in r] = 1 — Tq/Th as discussed in any basic 
textbooks. 



The validity of the second principle of thermodynamic 
reflects, as a corollary, the impossibility of realizing per- 
petual motion of the second kind. This undisputed ver- 
dict is recognized by the whole scientific community, de- 
spite the fact that the pursuit of perpetual motion re- 
mains a popular subject of speculations. If such a device 
existed, it could transform heat into useful work, with- 
out requiring any dissipative losses, a procedure that is 
ultimately against the second law. 

As already emphasized, the above picture is at core 
of our current understanding of physics, it appears to 
be solidly grounded on experiments and no doubt is 
cast on its correctness. Recently, however, a large body 
of evidence has emerged which testifies on the peculiar 
and rather unintuitive thermodynamic properties of long 
range interacting systems [j| . Microcanonical and canon- 
ical statistical ensemble can be incquivalent, an intriguing 
phenomenon that has been associated to e.g. negative 
(microcanonical) specific heat @, @] and susceptibility 
(for magnetic-like systems) [1, |]|. In classical thermo- 
dynamics, the specific heat measures the energy required 
to increase the temperature of a unit quantity of a sub- 
stance by a unit of (thermodynamic) temperature. As 
common sense dictates, providing energy to the system 
should induce heating, meaning an increase in the system 
temperature. At odd with intuition, non-Boltzmannian 
systems (as those displaying long range traits) can cool 
down when forced to increase their internal energy. The 
canonical specific heat is however always positive [Ioj |. 
a mathematical constraint that descends from statisti- 
cal mechanics. For non-Boltzmannian systems, however, 
the kinetic temperature differs from the thermodynamic 
one, the former providing a direct measure of microscopic 
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agitation. By elaborating on this scenario, it was demon- 
strated that negative kinetic specific heat can be 
measured in the canonical statistical ensemble. In fact, 
due to a complex self-organized interplay between col- 
lective and microscopic degrees of freedom, the particles 
belonging to the system under scrutiny can gain in kinetic 
energy when an energy quota is passed from the system 
to the heat reservoir. The opposite holds when the energy 
is instead released from the thermostat to the system. In 
reference fill ] . it was speculated that this peculiar pre- 
rogative traces back to the non-Boltzmanniann (equilib- 
rium) character of the particles' statistics. This could 
eventually open up the perspective of realizing novel, 
more efficient thermal devices. 

This paper is devoted to shed light onto such an is- 
sue, elucidating on a viable direction of investigation of 
paramount importance for both applied and fundamen- 
tal research. Working within a paradigmatic mean field 
model, and assuming the conventional definition of the 
canonical ensemble, we rigorously show that a thermo- 
dynamic device operating with a non-Boltzmannian fluid 
can be conceived, which violates the second law of ther- 
modynamics. The only possibility we see to reconcile 
our results with the second principle is to modify the 
Fourier law of heat conduction, by assuming that heat is 
transported along the gradient of the kinetic temperature 
and not the thermodynamic temperature, as classically 
believed. This modification is not innocent and can im- 
pact the upper boundary of thermal efficiency, as speci- 
fied by Carnot theorem. In both ways, we hence argue 
that working with unconventional non-Boltzmannian flu- 
ids can alter our current perception of the concept of effi- 
ciency, beyond the standard Boltzmann picture to which 
predictive tools are at present indirectly anchored. 

The long-range non Boltzmannian fluid. 

Let us now turn to introducing the general background 
that defines the playground arena for our thermody- 
namics derivation. We will deal with a simple one- 
dimensional model and label with 9 the configuration 
variable, p being its conjugate momentum. In the con- 
tinuum limit, we assume the single particle distribution 
function f(9,p,t) to obey to the following Vlasov equa- 
tion: 



2i ?L_wdi = 

dt P d9 dO dp 



(1) 



where the potential V(9)[f] is self-consistently defined in 
terms of the distribution itself. The problem is hence 
cast in a general setting which embraces diverse fields 
of application, from plasma physics to cosmology. Our 
declared aim is to construct and study an (equilibrium) 
cyclic machine working with a Vlasov fluid. To this end, 
we shall specialize further by assuming a specific form of 
the potential, which will enable us to carry out the sought 
calculations in details. The implications of our analysis 
can be however extended to other realms, provided the 
prototypical equation ([IJ holds true. In fact, equation 
([T]) yields a Fermi like entropy functional, which embeds 



dynamical constraints and that can be derived from first 
principles through the usual counting of accessible micro- 
scopic states. As a test case, we consider a paradigmatic 
example of long range systems, the so called Hamiltonian 
Mean Field (HMF) model, subject to an external mag- 
netic forcing of strength h. The HMF is a popular model 
of spins, interacting all-to-all via a cosinus like potential. 
V(9)[f] takes in this case the form: 



V{9)[f] = 1-M[f] cos(0) + hcos(9) 



(2) 



where M[f] stands for the magnetization: 

/7T POO 
/ f (9, p,t) cos 9d9dp , (3) 
-7T J — OO 

which measures the degree of spatial spins' bunching. In- 
cluding the kinetic contribution, the specific energy reads 
4.f] = I J(p 2 /2)f(9,p,t)d9dp - (M 2 - l)/2 + hM. In 
its original formulation, the HMF is a discrete model 
composed of N particles (spins) in mutual interaction, 
and has been thoroughly studied in the unperturbed 
limit h = 0. Starting from out-of-equilibrium initial 
conditions, the system gets trapped in Quasi-Stationary 
States (QSS) @, whose lifetime diverges when increas- 
ing the number of particles N. When performing the 
thermodynamic limit (TV — > oo) before the infinite time 
limit, the system remains permanently confined in the 
QSS, which are characterized by non Gaussian distri- 
butions and hence substantially differ from the deputed 
Boltzmann-Gibbs equilibrium. It has been shown that 
QSS's can be related to the stable (equilibrium) steady 
states of the Vlasov equation ([1]) describing the system in 
the continuum limit [f| [HI Studying the equilibria of the 
Vlasov equations could be hence equivalently interpreted 
as addressing the out of equilibrium QSS dynamics of the 
associated HMF model. 

The microcanonical (constant energy) behavior of sys- 
tem ^ can be investigated via a direct statistical 
mechanics treatment, inspired to the seminal work of 
Lynden-Bell [l2|, and originally developed in the con- 
text of stellar dynamics (see, e.g., [HI). Starting from an 
arbitrary initial condition, the Vlasov evolution makes 
the single-particle distribution function stir and distort 
in phase space. The hypervolumes of its values 77 (lev- 
els of phase density) are however conserved. Label with 
p(9, p, 77) the probability density of finding the level of 
phase density 77 in a small neighborhood of the position 
(0,p) in phase space. One can thus define a "coarse- 
grained" distribution function f(6,p) = / p(9,p,r))rjd7]. 
The microcanonical equilibrium distribution maximascs 
the mixing entropy s = — J p(9, p,ij) In p(9,p,rj)di']d9 dp, 
while conserving the invariant of the motion, as pre- 
scribed by standard statistical mechanics. For the sake 
of simplicity, and to allow for relatively straightforward 
analytical progress, we solely consider the case of a two- 
level water-bag initial conditions, i.e. r\ = /q or 77 = 



3 



|25[. Within this approximation, the entropy reads: 



s[f] 



dpde 
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(4) 

where / = fop. By maximizing the entropy functional, 
one obtains the microcanonical solution in terms of a 
fermionic distribution (see below), which coincides with 
a particular class of stable stationary solutions of the 
Vlasov equation. 

Given the microcanonical entropy (that self- 
consistcntly depends on the energy e), one can 
always calculate the canonical rescaled free energy, 
3, /), as the Legendre-Fenchel Transform (LFT): 
5, /) = inf e [fte — s(e, /)] . This is a universally 
adopted procedure: ft is a free parameter equal to the 
inverse of the (constant, in the canonical ensemble) 
thermodynamic temperature of the system. Requiring 
the free energy to be stationary, under the dynamical 
constraints, one recovers a closed mathematical pre- 
diction for the particles distribution, which again is 
fermionic. Concretely, the canonical distribution is: 



f(0,p) 



fo 



1 + e Pfo(p /2 - M[f] cos6-h cos 9) + a 



(5) 



where a plays the role of a Lagrange multiplier as- 
sociated to mass conservation. The latter equation is 
formally identical to the one obtained by requiring the 
entropy Q to be stationary in the framework of the mi- 
crocanonical ensemble [13] ■ In that case, ft and a are 
Lagrange multipliers associated to energy and mass con- 
servation. In both cases, the respective Lagrange mul- 
tipliers can be calculated by solving the self-consistent 
system of equations obtained by imposing the constraints 
conditions. Consequently, the coarse grained distribution 
function / and the magnetization M can be calculated, 
for any fixed initial condition, i.e. given /o, h and e (resp. 
ft) in the microcanonical (resp. canonical) setting. 

A remark is mandatory at this point, which will prove 
crucial in the forthcoming discussion. The thermody- 
namic temperature T t h is equal to 1/ft when working 
in the canonical framework, or specified by the funda- 
mental relation (<9s/<9e) _1 , if the microcanonical view- 
point is instead adopted. A system and a thermal reser- 
voir are said to be in equilibrium (no heat is exchanged) 
if their respective thermodynamic temperatures match. 
The canonical ensemble, as defined by the above LFT, 
provides hence the appropriate descriptive scenario for 
the system dynamics, ft being fixed as the inverse of 
the thermostat thermodynamic temperature. If a dif- 
ference in temperature manifests between two adjacent 
systems, the heat flows against the temperature gradient 
so to eventually restore the equilibrium condition. This 
is the celebrated Fourier's law. With reference to both 
microcanonical and canonical pictures, one can alterna- 
tively consider the definition of kinetic temperature [l5j , 



namely = J p 2 f(9,p)d8dp. For a Boltzmannian sys- 
tem f(0,p) tx exp(— ftp 2 1 2) [l6[ and it is therefore imme- 
diate to see that = l/ft = Tth- When instead the sys- 
tem displays non-Boltzmannian traits, Tk ^ T t h- Unintu- 
itive phenomena can thus set in as the emergence of neg- 
ative kinetic heat in the canonical ensemble, a possibility 
demonstrated in [ill ]. Moreover, for a non-Boltzmannian 
fluid, a word of caution should be exercised on the usage 
of Fourier's heat law: does the heat flow across the ki- 
netic or thermodynamic temperature gradient? We shall 
return later on this crucial point. 

Constructing the thermodynamic cycle. Having 
reviewed the definition of canonical and microcanonical 
ensembles, we can formally construct a thermodynamic 
cycle, reminiscent of Carnot's one, for an hypothetic de- 
vice working with the Vlasov fluid introduced above. We 
will in particular consider two isothermal transforma- 
tions, the reservoirs being characterized by the (thermo- 
dynamic) temperatures Th (hot) and Tq (cold). The 
cycle is then closed by assuming two adiabatic transfor- 
mations. As concerns the latter transformation, we recall 
fl7| that the first principle of thermodynamics reads: 



5Q = de- hdM. 



(6) 



Here, Q is the heat provided by the environment to the 
system. The adiabatic transformation SQ = is hence 
realized by tuning the field strength as h = de/dM 
along the path. A representative cycle obtained by im- 
plementing the procedure outlined above is depicted in 
figure [TJ Each point displayed in the three dimensional 
space (h, M, e) is either a stationary solution of the en- 
tropy functional s (microcanical, adiabatic lines) or of 
the free energy <j> (canonical, isothermal). Notice that 
one needs to circulate the cycle in figure [1] (lower panel) 
counter-clockwise, so to extract a positive work amount 
(Lhmf = — f cyc ie hdM), The heat exchanged with the 
hot and cold reservoirs, labeled respectively Qh and 
Qc, can be readily evaluated as an application of the 
first principle (JB]) and upon calculation of the work per- 
formed along the isothermal tracks. Surprisingly, and as 
a straightforward application of the above procedure, we 
obtain Qh < and Qc > 0. Moreover, Qc > \Qh\- This 
is at odds with the physical constraints imposed by the 
second law of thermodynamics. The virtual machine that 
we have constructed works in fact as a refrigerator, cool- 
ing down the reservoir with the lower temperature. In 
doing so, it however generates a positive work load, which 
violates Clausius formulation of the second law of ther- 
modynamics. In principle, and admitting this conclusion 
to be correct, one could insert between the two considered 
reservoirs a standard thermal machine that, in a cycle, 
receives from the hotter bath the heat quota Qh (he. the 
very same released by our Vlasov device) , and returns to 
the colder reservoir the heat Q' c , such that \Q' C \ < Qh, 
while producing a positive work L. The machine ob- 
tained by pairing together the two considered devices, 
subtracts heat from a single reservoir (Qc — \Q'c\ > 0) 
and returns the positive work L to t = Lhmf + L, against 
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the second law of thermodynamics in the Kelvin-Planck 
formulation. In other words, the machine that we have 
assembled, and which intimately exploits the unintuitive 
features of long range systems, is characterized by an ef- 
ficiency r] = L t ot/{Qc — \Q'c\) ec l ua l t° one - This is so 
unexpected that, at first sight, it would be tempting to re- 
ject the results as a whole. It is however possible to elab- 
orate on a thermodynamic interpretation of the results 
and pinpoint to the specific, and highly unconventional 
phenomenon, that is ultimately responsible for such an 
astonishing behavior. Then we will come back to dis- 
cussing the foundation of the theory, especially the cus- 
tomarily adopted definition of canonical ensemble, and 
suggest that the validity of the second principle could 
be eventually restored provided the kinetic temperature 
and not the thermodynamic one is employed in the law 
of Fourier for the heat transport. 

Negative kinetic heat capacity and the violation 
of the second principle of thermodyamics. 

The cycle reproduced in figure [1] extends within a do- 
main of the parameters space yielding magnetized states. 
In , it was shown that negative (kinetic) specific heat 
is found in the canonical ensemble, i.e. along the isother- 
mal lines of the cycle, when the magnetization is different 
from zero. This fact is indeed responsible for the appar- 
ently unphysical conclusions reached above. Consider in 
fact the first law and isolate the two terms that enter 
the right hand side of equation (|6]). The second term 
is related to the actual physical work: when proceeding 
counter-clockwise along the isothermal transformation, 
in contact with the hot reservoir, the quantity —hdM is 
positive as can be deduced by looking at figure [T] (lower 
panel, red line). In practice, and as it happens for an 
ideal gas performing an isothermal transformation, the 
work applied perturbes the system from its equilibrium 
condition, while constraining it to the same thermody- 
namic temperature of the bath. A (positive) heat quota 
{—hdM) is consequently transferred from the bath to the 
system, when circulating counter-clockwise. The oppo- 
site clearly holds when the the cold isothermal is con- 
sidered, together with the chosen versus of circulation: 
—hdM is negative, see figure [T] (upper panel, blue line). 
While the thermodynamic temperature is kept constant, 
the kinetic one changes along the considered transfor- 
mation: in particular, Tk increases along the isothermal 
at Th, decreases when the system is evolved in contact 
with the colder bath Tc- 

The other contribution to the heat SQ, the first term 
(de) on the right hand side of equation ©, gets also mod- 
ulated along the transformation. In figurc[2J the energy e 
is shown to decrease when Tk increases, for both cold and 
hot isothermals (resp. red and blue). The kinetic specific 
heat is thus negative, a condition that realises because of 
the non-Boltzmannian character of the system. Increas- 
ing the kinetic temperature, as the system does when in 
contact with the hot bath Th , implies decreasing the en- 
ergy e, and so contributing with a negative term to the 
global heat exchanged with the reservoir. Incidentally, 
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FIG. 1: Upper panel: The thermodynamics cycle is repre- 
sented in the three dimensional space (h, M, e). The red line 
stands for the isothermal transformation carried out at Th, 
while the blue line refers to the choice Tc < Th- Here 
T H = I/Ph = 0.33 and T c = l//3c = 0.26. The adia- 
batic transformations that close up the cycle are depicted 
in green. The heat exchanged with the reservoirs are re- 
spectively: Qh — —0.0363 (heat flows out the system) and 
Qc = 0.0416 (the heat flows towards the system). Here, 
/o = 0.12. Lower panel: The cycle is projected in the (M,h) 
subspace. The arrows indicate the direction of circulation 
that produces a positive work. 



the latter contribution turns out to be larger in absolute 
value, than the positive contribution associated to the 
magnetic work. The total contribution is hence negative 
and for this reason the system cedes heat to the hot reser- 
voir. Similar conclusions can be drawn when looking at 
the other isothermal transformation: now de is positive 
(the kinetic temperature decreases) and overcomes the 
negative —hdM contribution. When kinetic and ther- 
modynamics temperature are equal (as they should for 
conventional Boltzmannian short range systems) the (ki- 
netic) specific heat has to be necessarily positive, and the 
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FIG. 2: The energy e of the system is plotted versus the 
kinetic energy Tk, along the two isothermal transformations 
(in red Th = 0.33, in blue Tc — 0.26). The transformations 
proceed along the directions highlighted by the arrows. Here, 
fo = 0.12. 



above surprising condition cannot realize. 
Discussion. 

Summing up we have analytically shown that non- 
Boltzmannian systems can apparently violate the second 
law of thermodynamics, and explained such an aston- 
ishing ability in terms their peculiar statistical mechan- 
ics, i.e. existence of negative kinetic specific heat in the 
canonical ensemble. Importantly, and direct con- 
sequence of the above analysis, a cyclic thermal device 
working with a non-Boltzmannian fluid can be imagined 
which has efficiency equal to one. It is however hard to 
accept such a conclusion, which would revolutionize our 
current understanding of the laws of nature. 

A way out in our view exists which potentially allows 
one to respect the second principle and confirm its valid- 
ity for the class of systems being here considered. The 



critical point has to do with the definition of the canon- 
ical ensemble, that we have here mutuatcd from equi- 
librium classical statistical mechanics, by performing the 
Legendre-Fenchel Transform of the entropy functional. 
In doing so, one requires that the canonical dynamics 
preserves the thermodynamic temperature of the system 
setting it to a constant value imposed by the external 
reservoir. Alternatively, one could imagine that non- 
Boltzmannian systems, when placed in contact with a 
thermal bath, would organize so to freeze their kinetic 
temperature to the temperature 26] of the bath. As a 
natural consequence, already stressed before, the Fourier 
law [TH should require that the heat flows along the ki- 
netic temperature gradient, and not along the thermody- 
namic one, as the definition of canonical thermodynam- 
ical equilibrium would indirectly imply. Pushing further 
this proposal, we could possibly ends up with a cyclic 
non-Boltzmannian engine which meets the second law re- 
quirement, while still exceeding the Carnot upper bound 
of efficiency, since the thermodynamic temperature is re- 
lated to the kinetic one by a non linear functional relation 
[IH . This is an important issue that touches the founda- 
tion of modern physics and that deserves to be carefully 
addressed. 



To conclude, we would like to stress that our study is 
relevant for a quite large class of systems and for many 
practical physical situations. In fact, the Vlasov equation 
is often invoked as the reference model for several exper- 
imental systems, involving a large number of individual 
components, e.g. charged systems, plasmas [191 1. 
electron lasers [20( and magnetic dipolar systems 
particular, it can be shown that, under some hypothe- 
ses, the specific system here considered can be formally 
reduced, in the limit for h = 0, to that governing the dy- 
namics of collective atomic recoil lasers (CARL) [22|, [23| ■ 
This important observation could open up the perspec- 
tive for an experimental verification of our theoretical 
predictions. 
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